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Abstract 

We propose that, in SU(3) gauge theories with fundamental quarks, confinement can be 
inferred from spectral density of the Dirac operator. This stems from the proposition 
that its possible behaviors are exhausted by three distinct types (Fig.l). The monotonic 
cases are standard and entail confinement with valence chiral symmetry breaking (A) or 
the lack of both (C,C’). The bimodal ( anomalous ) option (B) was frequently regarded 
as an artifact (lattice or other) in previous studies, but we show for the first time 
that it persists in the continuum limit, and conclude that it informs of a non-confining 
phase with broken valence chiral symmetry. This generalization rests on the following. 

(а) We show that bimodality in Nj=0 theory past deconfinement temperature T c is 
stable with respect to removal of both infrared and ultraviolet cutoffs, indicating that 
anomalous phase is not an artifact. (/?) We demonstrate that transition to bimodality 
in N /=0 is simultaneous with the loss of confinement: anomalous phase occurs for 
T c <T < T c h, where T c h is the valence chiral restoration temperature. ( 7 ) Evidence 
is presented for thermal anomalous phase in Nj=2+1 QCD at physical quark masses, 
whose onset too coincides with the conventional “crossover T c ”. We conclude that the 
anomalous regime T C <T < T c is very likely a feature of nature’s strong interactions. 

(б) Our past studies of zero-temperature Nj=12 theories revealed that bimodality also 
arises via purely light-quark effects. As a result, we expect to encounter anomalous 
phase on generic paths to valence chiral restoration. We predict its existence also for 
Nf massless flavors (T = 0) in the range N)- < Nj- < N^ 1 , where N)• could be quite low. 
Conventional arguments would associate N c J l with the onset of conformal window. 

1. Introduction. SU(3) gauge theories with quarks in fundamental representation are 
important in elementary particle physics. Indeed, the strong dynamics of “real-world” quarks 
and gluons is believed to be described within this context, and the near-conformal infrared 
behavior in theories with many massless flavors (just below conformal window) may be 
relevant for theories of technicolor type. In addition, the systems within conformal window 
are attractive test beds of conformal dynamics in four dimensional quantum held theory. 

Since the quarks of nature are not massless, it is desirable to consider all theories of the 
above type, i.e. including any number of arbitrarily massive quarks. Denoting such theory 
space at zero temperature as T 0 and at any temperature as T D T 0 , this wide landscape 
involves various kinds of dynamics, expected to be well distinguished by corresponding vac- 
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Figure 1: Confinement/vSChSB structure of T via infrared behavior of Dirac spectral density. 

uum properties. [] In that regard, confinement and spontaneous chiral symmetry breaking 
(SChSB) dominate the thinking about strongly interacting dynamics. Nevertheless, classi¬ 
fying elements of 7 in these terms involves both conceptual and practical issues. 

One problem is that SChSB is well-defined only when at least a pair of quarks is massless. 
However, the vacuum of any theory in 7 can be probed for its ability to support long-range 
order via valence Goldstone pions (see e.g. BED- Formally, a pair of valence quarks and 
a pair of action-compensating pseudofermions is added to the system [3]. In the massless 
valence limit, flavored chiral rotations of valence fields rj, f) are elevated to symmetries, which 
the vacuum either respects or not. The latter possibility entails valence spontaneous chiral 
symmetry breaking (vSChSB), indicated by non-zero value of ( fjr /). SChSB and vSChSB are 
identical notions whenever SChSB is meaningful, but vSChSB provides for a symmetry-based 
distinction among vacua over the whole set 7. 

The situation is different in case of confinement, whose meaning is quite intuitive through¬ 
out T, but a satisfactory consensus on its precise definition in such general context (even 
within To) is lacking. Moreover, existing definitions are technically difficult to verify in a 
given theory. Convenient symmetry-based interpretation is only available when quarks do 
not affect vacuum correlations: in pure glue (N/=0) theories at finite temperature [41 E|. 
In this restricted case, the breakdown of Z 3 symmetry, signaled by non-zero Polyakov line 
(L) 7^ 0 or its magnitude (\L\) > 0, distinguishes confined and deconhned dynamics. 

Despite of issues with its formal definition over T, confinement can have simple and unam¬ 
biguous dynamical signatures. Indeed, since the phenomenon shapes dynamics in a crucial 
manner, it is reasonable to expect that even simple vacuum observables are qualitatively 
affected. While details of this influence depend on the unknown specifics of confinement 
mechanism, the logic here is that of reverse engineering: if simple signatures can be empiri¬ 
cally developed, consistently with existing results of lattice simulations and with theoretically 
clean N/=0 case, then a valuable insight into the mechanism is obtained. In addition, prac¬ 
tical benefits for investigating the phase structure of T are likely to be acquired. 

Following such rationale, we propose that vSChSB/conhnement structure of T may be 
inferred from low-energy behavior of Dirac spectral density p(A) in these theories. Gauge 
invariant p(A) is thus viewed as a vacuum object already known to encode vSChSB via its 
“infinitely infrared” behavior (p(A—>-0) in infinite volume) and Banks-Casher equivalence [6j. 
The suggestion is that properties of p(A) in finite infrared regime reflect vSChSB/conhnement 
combination as shown in FigjlJ case (A) corresponds to vSChSB with confinement, (B) to 

1 Note that, in what follows, we will not differentiate between “vacuum” (zero-temperature concept) 
and “equilibrium state” (finite-temperature concept). Path integral description of the state (via typical 
configurations) makes little difference between the two at the technical level. 
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vSChSB without confinement, and (C,C’) to valence chiral symmetry and no confinement. 
Note that option (A) represents generic behavior observed in low-temperature “real-world” 
QCD simulations, and (C,C’) are standard possibilities for chirally symmetric vacuum]^] We 
refer to the atypical case (B) as anomalous. There are few explanatory comments to make. 

(i) Dependencies of Fig jl] refer to limits approached in asymptotically large volumes. 
Non-monotonicities arising e.g. in small volumes at fixed topological charge are not invoked. 

(li) While Banks-Casher relation is a kinematical constraint, the above connection be¬ 
tween vacuum properties and p(A) is dynamical, arising due to the nature of strong force. 

(in) Note that confinement without vSChSB is not among the possibilities in T, aligning 
with conclusions of Ref. [7j. 

(tv) We implicitly assume a fully chiral lattice definition of valence quark dynamics via 
overlap Dirac operator [8]. This makes vSChSB and the proposed classification of phases 
sharply defined at the regularized level. Note that the nature of phase identification is such 
that it is immaterial whether bare or renormalized [9] spectral density is used. It is equally 
immaterial that p(X — y 0 ) (valence chiral condensate) may diverge even at the regularized 
level in some theories [TO] : this just means that the system is necessarily in phase (B). 

(v) Bimodality in p(A) was first observed on N/=0 backgrounds [TTj above deconfinement 
temperature T c . However, its cutoff dependencies have not been systematically studied 
amidst suspicions (see e.g. [12]) that it may be a lattice artifact. Consequently, substantiating 
the existence of anomalous phase in 7 is central to our discussion. 

(vi) Recent observation of anomalous p(A) in zero-temperature dynamics with many light 
flavors mm provided an important impetus for current generalization: since both thermal 
and light-quark effects induce anomalous phase, it is reasonable to expect its occurrence 
along generic paths in 7 leading to valence chiral symmetry restoration. 

(vii) In the current context, the only role of light valence quarks is to probe the vacuum: 
while not influencing it in any way, their infrared dynamics is sensitive to its long-range 
properties. Contrasted with type (A), the anomalous density (B) then signifies a different 
mechanism for Goldstone-like correlations, relevant when confinement is absent. 

In what follows, we first discuss new results from lattice QCD needed for this proposal, 

1. e. points (a)-( 7 ) outlined in the Abstract. We then proceed to generalize (Fig. [ 6 ]), taking 
into account point (h), and suggest the existence of anomalous phase preceding the conformal 
window. Selected aspects of our results, conclusions and their consequences, are elaborated 
upon in the last section. 

2. Reality of the Anomalous Phase. Basic requirement for viability of the proposed gen¬ 
eral scenario is establishing the existence of deconfined dynamics with vSChSB and anoma¬ 
lous spectral density anywhere in 7. Thermal N/=0 theory is ideal for this purpose since 
all three elements involved are well-defined and lattice-accessible. Working with Wilson’s 
lattice theory at T/T c =1.12 (fixed by r$T c [H]), where we observed bimodality in spectral 
density previously EE!, our aim is to perform crucial stability tests with respect to infrared 
and ultraviolet cutoffs. Our preliminary work in this direction was reported in Ref. pRJ]. The 
results below refer to Z 3 -broken vacuum with “real Polyakov line” ra. which gets selected 
when deforming the theory away from N/=0. 

2 Case (C’) involving strict “gap” in density is very difficult to infer from any practical simulation and 
may or may not exist in T. However, the distinction between (C) and (C’) is immaterial for our purposes. 
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Figure 2: Anomalous phase in Nj=0 theory at T = 1.12T C and lattice spacing a = 0.085fm. 


Overlap Dirac operator (p = 26/19) constructed from its Wilson counterpart (r = 1) is 
used to define valence quark dynamics throughout this work. Spectral density p(A) is the 
right derivative of cumulative function cr(A) = ( ^ 0<Ai<A 1 )/V, where V is the 4-volume and 
A i (real numbers) have magnitudes of Dirac eigenvalues and signs of their imaginary parts. 
Indicated exclusion of exact zeromodes from counting is harmless since their effect vanishes 
in the infinite-volume limit. Given the finite statistics of any simulation, coarse-graining 
of p(A) is unavoidable. In the absence of suspicion for singularity, we use the symmetric 
definition away from origin, namely 


p(M) s 


A > <5/2 


( 1 ) 


When estimating the infinite-volume value of p(A —> 0) from finite systems, it is convenient 
to work with “right derivative” form, i.e. p(A = 0, A) = cr(A)/A. Valence chiral condensate 
with overlap is identical to tt lim A ^. 0 limy-voo p(A, V), as is formally in the continuum. 

Focusing first on the infrared, we simulated N 3 x 7 systems at N = 20, 24, 32,48. To set 
T = 1.12T C requires lattice spacing a = 0.085fm (/5 = 6.054, ro = 0.5fin), giving linear size 
L — 4.1 fm to the largest system. Low-lying overlap eigensystems were computed and bimodal 
p(A), such as in the top left plot of Fig. [2j was found at all volumes. The spectrum exhibits 
excellent volume scaling, exemplified in the bottom left plotj^J except for the most infrared 
point, which follows a growing trend ensuring the anomalous shape (B) in the infinite-volume 
limit. The coarse-graining 5 — 20 MeV is sufficiently fine except for the most infrared bin 
(0,20)MeV. Here p(A) changes rapidly, as the closeup with 5 = 4 MeV (top middle) shows. 
Further lowering of 5 (top right) only affects the most infrared point, and we have thus 
arrived at the resolution revealing the shape of the anomalous peak. Importantly, this shape 
is also stable under the change of infrared cutoff: comparing N = 32 to N = 48 (lower middle 
vs top middle), only accumulation in the lowest bin changes, growing again. We conclude 

3 The horizontal lines guiding the eye are fits to a constant from three largest systems. 
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Figure 3: Ultraviolet cutoff and anomalous behavior: N/=0 theory at T = 1.12T C , L = 2.0fm. 


that, in the infinite volume, p( A) has a positive (possibly infinite) local maximum at A = 0, 
following behavior (B), with anomalous peak of finite width (few MeV). 

While anomalous behavior (B) alone implies it, we address valence chiral symmetry 
breaking explicitly (lower right plot of Fig. [2j) via A-dependence of p(A = 0, A). Note that, 
in any finite volume, liniA^o p(0, A, V) = 0, and the associated downturn, shown for N = 32, 
will occur at sufficiently small A. However, as expected in broken theory, the break is pushed 
toward zero with increasing volume. Moreover, the estimates at any fixed A are growing, 
making vSChSB all but certain. At T — 1.12 T c the system is deconfined, and thus all elements 
of the proposed connection are in place for this case. 

Turning to ultraviolet cutoff, we fixed the physical volume of A = 24 theory (L = 2.04fm), 
and drove the system toward continuum limit, while keeping T/T c = 1.12. This resulted 
in sequence N x N t = 24 x 7, 28 x 8, 34 x 10, 42 x 12 at a = 0.0850,0.0744, 0.0595, 0.0496 fm 
respectively. The choice of moderate physical volume was chiefly motivated by computational 
considerations. While data presented in Fig. [2] indicate volume effects in the infrared, the 
qualitative behavior is still very clearly anomalous at this volume, and the setup allows us 
to explore the system rather deep into the continuum limit which is essential. 

The bimodal p(A) was indeed found at all cutoffs, with global view for the system closest 
to the continuum displayed in Fig. [3] (top left). The associated closeup with proper resolution 
is also shown (top middle). Note that, contrary to changing the volume, varying ultraviolet 
cutoff is not expected to affect very infrared scales significantly, unless a qualitative change 
of dynamics (phase transition) occurs, separating lattice and continuum-like behaviors. It 
is decisively the former scenario that is observed, and illustrated by comparisons to the 
situation at smaller cutoff (lower left and middle plots). The scaling of first and second data 
point at <5 = 20 MeV (top right) indicates that the bimodal shape of p(A) will be preserved in 
the continuum limit at this resolution. Moreover, the accumulation of very near-zeromodes, 
i.e. p(0, A) with resolution well within the natural width of the anomalous peak (A = 4 MeV), 
is stable with cutoff (lower right plot), offering no hint of a qualitative change. 
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Figure 4: Thermal transitions to deconfined and anomalous phase coincide in N/=0 theory. 


Given the demonstrated dramatic volume effects leading to anomalous phase (B) at 
physically relevant cutoff, and the stability of bimodality at fixed volume under ultraviolet 
cutoff, we conclude that continuum N/=0 theory at T/T c = 1.12 is in the anomalous phase. 

3. Confinement and Anomalous Spectral Density. An important aspect of the pro¬ 
posed general scheme is that the transition from confined to deconfined theory in T coincides 
with the transition from regular type (A) of p(A) to anomalous type (B). To show this for 
N/=0 thermal transition, we selected a system just below T c . Replacing N t — 7 of our volume 
study with N t = 8 (keeping (3 = 6.054, a = 0.085fm) corresponds to temperature T = 0.98 T c . 
This nominal assignment uses universal continuum value r 0 T c = 0.7498(50) of Ref. [2], and 
T/T c at finite cutoff is expected to shift slightly upwards. We thus checked the behavior 
of Polyakov loop in large volumes up to N = 48, and confirmed that the system at hand is 
indeed confined and at the very edge of the transition. 

Rg.g (left) shows spectral density for theory so tuned, on 20 3 x 8 lattice. Type (A) 
behavior is found, as predicted, with flat dependence at low energies indicating a system 
on the verge of the transition. To compare this to the situation at temperature just above 
deconfinement point, we increased the coupling slightly (/3 = 6.0783, a = 0.0817 fm) to set 
T = 1.02T C . Spectral density for such theory on 36 3 x 8 lattice (middle) shows a dramatic 
change to sharp bimodal behavior: there is little doubt that the Z 3 deconfinement transition 
coincides with the transition to anomalous phase. Note that we also show the result for 
system at T = 1.12T C (32 3 x 7, /3 = 6.054) in a volume comparable to T = 1.02T C case, 
exemplifying the situation well inside the anomalous region (right). 

4. Thermal Anomalous Phase of “Real World” QCD. The above evidence of anoma¬ 
lous phase in Nj=0 QCD, and its association with deconfinement, gives an essential impetus 
for its existence in other corners of T. It is of particular interest to clarify whether anoma¬ 
lous dynamics is part of physical reality for nature’s quarks and gluons. We address this 
via N/=2+l QCD at physical quark masses, i.e. light mass of (rn u + m^)/2 and heavy of 
m s , since it is well established that this theory provides for very precise representation of 
relevant strongly interacting physics. Among other things, calculations in this context led 
to the conclusion that thermal transition of strong interactions is a crossover da 

To probe the crossover region, we utilize N t =8 lattice ensembles of Wuppertal-Budapest 
group (see [18] and references therein), used in precise determination of transition temper¬ 
atures. On a technical side, the simulated theory involves Symanzik-improved (tree level) 
gauge action and stout-improved staggered fermions. The physical line of constant physics 
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Figure 5: Anomalous phase across the thermal crossover in N/=2+l QCD at physical point. 

was defined via fixing /a/ 777 -^ Ik/^k and both scaling violations and taste-splitting stag¬ 
gered fermion effects were shown to be small at gauge couplings involved. Due to its crossover 
nature, transition temperature is not a unique concept and depends both on the observable 
used and the defining condition chosen. For reference, the temperatures associated with 
inflection points of (light) scalar density (“condensate”) and Polyakov line were reported as 
T c fa 155 MeV and 170 MeV respectively. Given that, we selected ensembles at temperatures 
T = 150, 175 and 200 MeV, to be close to both ends of the crossover, and to examine the 
possible presence of anomalous phase past the transition region. 

The resulting overlap spectra are shown in Fig. [5] On the lower edge of the crossover 
(left plot), spectral density is extremely flat in the infrared, but follows the standard (A) 
behavior. On the other hand, just past the Polyakov line crossover temperature (middle 
plot), the anomalous behavior clearly sets in. There is thus little doubt that the appearance 
of anomalous phase closely follows the conventional measures for expected transition to 
deconfinement, i.e. conventional “T c ”. At T = 200 MeV, well above the established transition 
region, the anomalous phase fully develops, featuring prominent anomalous peak and the 
large degree of depletion at intermediate scales. 

Given the extensive checks performed on the ensembles used, both with respect to contin¬ 
uum limit and sufficiency of volume^] (see [IB] and references therein), we conclude that the 
anomalous phase T C <T < T c h, i-e. deconfined phase with broken valence chiral symmetry, 
exists in high-temperature dynamics of strong interactions, and follows the general scheme 
proposed here. 

5. Generic Anomalous Phase, N£ and NjP. We have presented evidence that anoma¬ 
lous phase with claimed properties exists in “T, as well as evidence that thermal effects lead 
to bimodal p( A) rather generically. Regarding the latter, apart from results presented here, 
observation of bimodal behavior in SU(N) theories was reported under varied but always 
thermal circumstances, e.g. in Refs. mmmmmmmm- Important additional 
ingredient is provided by the recent finding that presence of light dynamical quarks, without 
any thermal agitation, also leads to anomalous phases at sufficiently large number of fla¬ 
vors mm- Thermal and light-quark effects thus appear to be analogous in this regard. The 
relevance of this stems from the fact that thermal effects (increasing T), light-quark effects 
(decreasing masses, increasing Nf ), and their combinations, are the only available freedoms 
in T to facilitate the transition from broken valence chiral symmetry to its restoration. This 

4 The smallest linear size involved here is L=3.9 fm at T = 200 MeV. 
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Figure 6: The structure of set 7 proposed here (left) and the conventional one (right). 


leads us to propose that anomalous phases commonly occur on paths to chirally symmetric 
vacua. Our evidence on simultaneity of transitions to deconfinement and anomalous p( A) 
(in known cases) completes the rational basis for the proposed picture. 

Note that for N/=12 theories, the above scenario was seen [2, 13] to imply the existence 
of a mass m c below which the anomalous phase appears. The phase may either extend down 
to rn — 0 or to non-zero m c h < m c , namely the possible point of valence chiral restoration. 
Here we wish to point out another special case of the above general argument, concerning 
the “path” in 7 parametrized by the number of massless quark species N f at T = 0. Indeed, 
it is widely believed that the N f =2 case represents a confining, chirally broken theory with 
type (A) spectral density. At the same time, since the work of Banks and Zaks [2_2J, it is 
expected that there is a critical number of flavors N J < 16.5, such that theories in the window 
N J < N f < 16.5 are both asymptotically free, and controlled at low energy by a conformal 
infrared fixed point j^] Theories in the conformal window cannot generate low-energy scales 
(at least not below the scale of conformality), and the standard scenario is that N J marks 
the common transition to deconhned, chirally symmetric phase with p( A) of type (C). 

However, based on the proposed picture, we predict that when increasing the number of 
flavors from Nj=2, quark-gluon dynamics will first lose confinement (at N^), and only then 
chiral symmetry breaking (at N/), generating the anomalous phase for 

2 < N} < N/ < Nf = N f (2) 

The last equality holds if the current view associating the onset of conformal window with 
chiral symmetry restoration is valid, ft is important in this regard that N/=12 may be in 
the anomalous phase mm- Indeed, the setup with yet fewer flavors is even more conducive 
to the possibility of anomalous vSChSB surviving masslessness (SChSB), since less of a 
condensate-destructive light-quark effect is generated. could thus be quite low. 

Schematic structure of the set T, we are proposing, is shown in Fig. [6] (left). At T = 0 
(set To) there are theories with all three types of p( A) and associated vacuum properties. We 
refer to type (A) as QCD-like for this purpose. Type (C) is expected to contain theories from 
conformal window and e.g. theories without asymptotic freedom. Heating up a QCD-like 
system generically involves both deconfinement (T c ) and separately valence chiral restoration 


5 Note that we consider N/ to be an integer, while boundary values such as N y, N^ and N/ half-integers. 













( T c h ), while anomalous type (B) theory only undergoes the latter, and type (C) neither. 
Moving within 7 at fixed temperature involves analogous transitions. Note that other kinds 
of changes, unrelated to vSChSB and confinement, may occur within the three regions. Also, 
we do not imply that phases (A) and (C) are strictly disconnected, with anomalous phase 
occurring along every path between them: the intention is to convey a typical situation. 

6. Important Points. In this work, we have proposed a specific association between in¬ 
frared behavior of Dirac spectral density, and confinement/valence chiral symmetry breaking 
structure of SU(3) gauge theories with fundamental quarks (Fig. [Tj) . There are few points 
we wish to discuss and highlight in this regard. 

I. A significant conceptual step implied by our analysis and represented by Fig. [6j is 
that valence chiral symmetry breaking without confinement is a commonplace for SU(3) 
gauge theories with fundamental quarks: the schematic view on the left, with many new 
transitions, is proposed to be a good representation of the situation, as opposed to the 
standard view on the right. The existence of anomalous phase in T, shown here, is pivotal 
for this transformation. 

II. The proposed association of Dirac spectral density with vacuum properties turns it 
into a valuable tool for investigating the phase structure of 7. Indeed, recall that standard 
methods, such as those used in thermal lattice QCD, are “transitional” in nature, requiring 
the study of changes in the theory space. On the other hand, utilizing p(A) is as convenient 
as order parameter at hand: it provides a definite answer for any standalone theory. A 
single-number indicator of anomalous phase (“order parameter”) can be defined as follows. 
Let A an be the largest scale such that p(A, V —*oo) strictly decreases on [0, A an ). Then 

Oan = er(A an ) > o (3) 

i.e. a non-zero volume density of anomalous modes, identifies the theory in anomalous 
phasej^] This definition makes Q a n analogous in nature to Q, the volume density of chirally 
polarized modes, which was proposed to be a valid order parameter of vSChSB in Ei- 

III. The importance of concluded anomalous (deconfined with broken valence chiral sym¬ 
metry) phase in physics of quarks and gluons is hard to overstate. Indeed, the temperature 
regime involved is relevant for strongly interacting dynamics associated with the creation 
of plasma-like state currently studied at RHIC and LHC. While the bimodal aspect of p(A) 
didn’t remain completely unnoticed with regard to quark-gluon plasma phenomenology [23], 
the reality of anomalous phase, and the associated details, provide strong rationale for de¬ 
tailed exploration of compatible models. We will focus on such issues in a separate study. 

The emerging evidence for strong interactions generating two “high-temperature” phases, 
rather than one, complicates the associated physics in an intriguing way. This complexity 
can be seen in chiral polarization properties of Dirac modes. Indeed, our results imply that, 
past the thermal crossover, quark dynamics generates rich multilayered chiral structure in 
Dirac spectrum, shown in Fig. 17 of Ref. |2j (left plot). 

IV. If the proposed anomalous phase separating QCD-like behavior and infrared confor¬ 
mality materializes, it is expected to be instrumental for understanding of massless quark 
dynamics, and of transition to conformality in particular. It may also complicate some op¬ 
tions being considered for realization of the walking technicolor scenario [2SE5]. While it 

6 Note that \ an itself is such an indicator. 
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is exactly the near-conformal theories that are of interest for that purpose, one should em¬ 
phasize that dynamics in anomalous phase is still strongly coupled, and that deconfinement 
doesn’t imply a complete absence of hadronic bound states. The theories from “anomalous 
window” could thus still be relevant in the context of walking technicolor. 

Note Added: Few days before the first version of this manuscript was released, the work 
relevant for point ( 7 ) of the Abstract also appeared [26]. While the discussion of [26] is 
carried out in the context of Ua{ 1) restoration problem, their results in N/= 2 +l theories 
with light mass not far from physical, are consistent with and corroborate our conclusion 
that thermal effects generically lead to the anomalous phase (as defined here). 
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